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Abstract 

We study the density of states measure for some class of random unitary band matrices 
and prove a Thouless formula relating it to the associated Lyapunov exponent. This class 
of random matrices appears in the study of the dynamical stability of certain quantum 
systems and can also be considered as a unitary version of the Anderson model. We 
further determine the support of the density of states measure and provide a condition 
ensuring it possesses an analytic density. 

1 Introduction 

The stability of quantum dynamical systems generated by time periodic hamiltonians is 
sometimes characterized by means of the spectral properties of the corresponding unitary 
evolution operator over a period, also called monodromy operator, see |Bel IHoll ICo3j . 
Unfortunately, even for this relatively simple time-dependence, except for certain specific 
models, e.g. |Ho2l iDFl IBo] . it is rarely the case that one has enough information about 
the actual monodromy operator so that a complete spectral analysis can be performed. 
Therefore, one resorts to different approximation techniques in some specific regimes to say 
something about the spectrum. For example, KAM inspired techniques, see e.g. |Be| ICol| 
ID^IADHIdUTv1IHY| . or adiabatic related approaches, see e.g. |B^IB^IHo4MN1| IJ1IN2|. 
have been used to tackle this problem. 

In case the complexity of the monodromy operator is important enough to forbid of a 
complete description of it, one may resort to a statistical modelization. It is the case in 
particular in the study of the quantum dynamics of electrons confined to a ring threaded 
by a time dependent magnetic flux, see e.g. the paper BBj and references therein. A mod- 
elization of this dynamics by means of an effective random monodromy operator taking 
into account the details of the metallic structure of the ring is considered and tested nu- 
merically in BB . Wc refer the reader to this paper and |BH.T| for a more detailed account 
of the construction of the monodromy operator. 

Motivated by this approach, the spectral analysis of a class of random and deterministic 
unitary operators, which contains the above monodromy operator, is performed in |BH Jj . 
The main characteristics of these unitaries is that, when expressed as matrices in some basis, 
they display a band structure: more precisely they are five-diagonal. While the construction 
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of the set of unitaries studied in |BHJj is patterned after the above mentionned physical 
model, we believe it may be useful for a wider class of problems. Moreover, in the regime 
we consider here, this set of unitaries bears strong resemblances with the Jacobi matrices 
related to self-adjoint discrete one-dimensional Schrodinger operators. 

Another motivation in that direction stems from the recent paper |CMVj where certain 
unitary infinite matrices associated to the construction of orthonormal polynomials on the 
unit circle are shown to display the same five-diagonal structure as our set of monodromy 
operators. These matrices are shown in |CMVj to be unitarily equivalent to unitary op- 
erators introduced almost ten years ago in |GT| for the study of a related trigonometric 
moment problem. Moreover, in the latter paper, some effects of randomness in the coeffi- 
cients of these operators were investigated. 

The goal of the present paper is to pursue the analysis of such random unitaries in 
the setting considered in the paper [BH Jj . The phases of the matrix elements of the five- 
diagonal operators are random variables and the deterministic modulus depend on one 
parameter only: if the phases are all set to zero, what we will call the "free case", the 
unitary operator depends on a "reflexion" coefficient r g]0, 1[, see below. However, while 
the analysis of |BHJj focused on spectral issues, i.e. proving singularity of the almost sure 
spectrum by means of a unitary version of the Ishii-Pastur theorem and the positivity of the 
Lyapunov exponent obtained via Furstenberg's Theorem, the main object of the present 
study of the density of states measure and its links with the corresponding Lyapunov 
exponent. 

More precisely, expressing the density of states as the density of eigenvalues of a series 
of unitary operators restricted to "boxes", we are able to state this relation as what is 
known as a Thouless formula. This formula allows to compute the Lyapunov exponent 
by means of the density of states and to recover the a.c. component of the density of 
states measure by means of a derivative of the Lyapunov exponent. A consequence of 
our version of Thouless formula is the extension of some results of |BHJ| providing, in 
particular, an explicit value of the Lyapunov exponent in these cases. We also prove the 
validity of the Thouless formula for the deterministic free case, by explicit computations of 
the relevant quantities. Taking advantage of the analogy of our unitary matrices with the 
one dimensional discrete Schrodinger operator, we characterize the support of the density 
of states in term of that of the distribution of the i.i.d. random phases. Finally, we provide 
an effective criterion ensuring analyticity of the integrated density of states in terms of the 
exponential decay rate of the Fourier coefficients of the the distribution of the phases. This 
result relies on some kind of propagation estimates for the free evolution. 

We note here also that a Thouless formula is proven for the unitary random operator 
studied by Geronimo and Teplyaev in |GT| . The corresponding random matrix V u is defined 
in the canonical basis of l 2 (1) as well but displays a different structure: for all k € Z, the 
vector V^ifk has non zero coefficients along ipj, for j = — do, • • • , k+1 only. Such matrices are 
also called Hessenberg matrices. However, the operator under consideration here presents 
different characteristics from the one of |GT| . or |GMV| . in particular regarding the way 
randomness appears in the phases of the matrix elements. 

The plan of the paper is as follows. Section 2 is devoted to the definition of the model 
and its basic properties. The density of states is introduced in the next section and Thouless 
formula is proven in Section 4. The statements about the support of the density of state 
and ist analyticity properties are made in Section 5, whereas an Appendix contains some 
technical items. 
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2 The Model 



We present here the unitary matrices we will be concerned with and recall some of its basic 
properties to be used later. 

The unitary operator we consider has the following explicit form in the canonical basis 
Wk}kez of l 2 {1) 

U u <p*k = irte^k(p 2k -. 1 + r*e- ir &><pz k 

+ irte'^+i ip 2k+1 - iV ir? 2fc+i ip 2k+2 

U u <pzk+l = -^e-^^k-i + itre-^^ 

+ rV^+i^fc+i +irte~~ ir >M+i(p 2k+2i (2.1) 

for any k G Z. According to |BH Jj . the random phases {f] k } k ei, are functions of some 
physically relevant i.i.d. random variables {(0 k , a k )} k& % on the torus given by 

rK = ffg + et-i + c%-c%_ 1 , (2.2) 

for all k E Z and the coefficients r,t g]0, 1[ are interpreted as reflexion and transition 
coefficients linked by r 2 + t 2 . We will identify the operator and its matrix representation 
(|2.1|) . Let us recall that these parameters are assumed to be different from their extreme 
values and 1, because in case r = 1 t = the operator U w is diagonal and if 

r = <^=^ £ = 1, it is unitarily equivalent to the direct sum of two shifts. Let us finally 
mention that U u is constructed in section 2 of |BHJj as a product of two unitaries given 
by infinite direct sums of 2 x 2 unitary blocks. 

2.1 Ergodic properties 

More precisely, let us introduce a probabilistic space (O, P), where O is identified with 
{T z }, T being the torus, and P = ®k&iFki where F 2k = P and P 2 fc+i = Pi for any k G Z 
are probability distributions on T and T the u-algebra generated by the cylinders. We 
introduce the set of random vectors on (Q, J 7 , P) given by 

Pk = {9k,a k ) : n^T 2 , keZ, 

e% = uj 2k , a% = oj 2 k+i- (2.3) 

The random vectors {Pk}kez are thus i.i.d on T 2 . 

We denote by U u the random unitary operator corresponding to the random infinite 
matrix (|2.1|) . In analogy with Jacobi matrices describing the discrete Schrodinger equation, 
we will also denote the vector ip k by the site k, k E Z. 

Introducing the shift operator S on f2 by 

S(u;) k = cj k+2 , keZ, (2.4) 

we get an ergodic set {S J }j<=z of translations. With the unitary operator Vj defined on the 
canonical basis of Z 2 (Z) by 

Vjtfk = <Pk-2j,Vk G Z, (2.5) 

we observe that for any j G Z 

U S i u = VjU u Vf. (2.6) 
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Therefore, our random operator is a an ergodic unitary operator. Now, general argu- 
ments on the properties of the spectral resolution of ergodic operators E W (A), where A 
is a Borel set of the torus T, ensure that this projector is weakly measurable, as well as 
E®(A) = P^E LJ (A), where x = p.p., a.c. and s.c, denote the pure point, absolutely con- 
tinuous and singular continuous components, see |CL| . chapter V. The analysis performed 
in |BH.7| for the case where {(9%, a%)}k£Z are uniformly distributed on the torus shows 
that the a.c. component of the spectrum of U w is almost surely empty. 



2.2 Lyapunov Exponent 

Let us proceed by recalling some facts concerning the Lyapunov exponent. It is shown in 
BB and BHJ that generalized eigenvectors defined by 

= e iX ^, 

tp = ^Cktfk, c fc G C, A E C (2.7) 



in our unitary setting can be computed by means of 2 x 2 transfer matrices due to the 
structure of the matrix U w . They are such that for all k E Z, ( |BHJj ) 

^ \ =T(k)( C2k -A (2.8) 

,C2fc+l / \ C 2fc-l / 

where the randomness lies in the phases %(A) = rf£(X) defined by 

%(A) = Vk + A, (2.9) 

and 

T(k) u = - e -^-i(A) (2.10) 
T(k) 12 = i r - ^r*o*k-iW - l) 

T{k) 22 = ™ e i7l2kix) + ^ ^ e *'(nfcW-u»_i(A)) + 1 _ g-iranCA) j _ 
Note the properties 

T(k) = T(r, 2k (X), V2k ^(X)) (2.11) 

whereas detT(fc) = ^Wzk-mk-i) j s independent of A. 

Therefore, knowing e.g. the coefficients (cq,ci), we compute for any k E N, 

C2k ) = T{k)...T{2)T{l)( C A^mh 
c 2 k+i J \ ci / V c i 



*J = T(-k + l)- 1 ---T(-l)- 1 T(0)- 1 (^j = (2.12) 



C- 2 k 
. C-2A:^ 

The dynamical system at hand being ergodic and the determinant of the transfer matrices 
being of modulus one, we get the existence of a deterministic Lyapunov exponent j(e lX ), 
for any A E C, such that 

lim ±-\n\\<S>(k)\\= 7 (e lX ) a.s.. (2.13) 

fc^ioo \k\ 

Writing e iX = z E C\{0}, we also know from classical arguments, see e.g. (T'KSl. that 7 
is a subharmonic function of z. 
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3 Density of States 



Following the standard approach in the self-adjoint case, we start by a definition of the 
density of states by averaging over the phases and invoking the Riesz-Markov theorem. 
Then we relate the density of state with alternative definitions in terms of the density of 
eigenvalues of truncations of the original operator to l 2 ([M, N]), as N — M — > oo. 

Definition: The density of states is the (non-random) measure dk on T defined by 

[ f(e tX )dk(X) := E[<¥>o|/(£4>o) + (H/(^)¥>l)]A (3.1) 

for any continuous function f : S 1 — > C. 

The average over the (fo and ip\ matrix elements is motivated by the forms of the matrix 
(|2.1[) and shift (|2.4[) . Note also that this definition makes dk a probability measure. 

Now we turn to the definition of appropriate finite size unitary matrices constructed 
from (|2.1j) . There are several possible constructions suited to our purpose. Those we 
use below result from considering U w provided with boundary conditions at certain sites 
forbidding transitions through these sites. Although such an interpretation is not needed 
in the sequel, let us mention it can be seen in |BH.T| . There, a more general unitary matrix 
than (|2.1[) is considered, whose reflection and transition coefficients (r^ifc) may depend on 
the index k, whereas (j2.1jl is a special case with r^ = r and t^ = t. Imposing tjy = there, 
one gets that the matrix takes a block structure which decouples the sites with indices 
smaller than iV from those with indices larger than N. 

Let us drop temporarily the sub- and super-scripts to in the notation. Fix N £ Z 
and consider the unitary operator U 2N on / 2 (Z) obtained from the original operator U by 
imposing the following boundary conditions at the sites 2N. Let U 2N be defined by 1)2 
for k g {2N, 2N + 1} where 

V2N-1 = V2N = V2N+1 = V2N+2 = (3.2) 

and, for k £ {2N, 2N + 1} 

U 2N (f2N = it<P2N-l + rip 2 N 

U 2N ip 2 N+l = T27V+1 + itip 2 N+2- (3.3) 

Similarly, a boundary condition imposed at site 2N + 1 defines JJ 2N+l by (|2,lj) for k ^ 
{2N, 2N + 1,2N + 2, 2N + 3} where 

V2N+1 = 1]2N+2 = (3.4) 

and, for k & {2N, 2N + 1,2N + 2, 2N + 3} 

U 2N+1 V2N = irte~ iri2N ^2N-i + r 2 e- iri2N if2N + it^N+i 

U 2N+1 <p2N+i = -t 2 e- l ^^2N-x + irte- ir >™<p2 N + r<p 2N+1 

U 2N+1 ^ 2N+ 2 = r<f2N+2 + irte-^+z^N+s - t 2 e~^+^2N+A 

U 2N+ \ 2 N+z = +it<f2N+2 + r 2 e-^^+V27V+3 + irte-^+^2N+4. (3.5) 

For any M G Z, the corresponding operator U M has a the block structure mentionned 
above and it is unitary. Then, given (M, N) G Z 2 such that M + 4 < N, one defines a 
unitary matrix JJ M,N on Z 2 (Z) by imposing boundary conditions at sites M and N. By 
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construction, U M ' N contains an isolated (N — M) x (N — M) unitary block on l 2 ([M + l, N]) 
we denote by V M ' N . Introducing the characteristic function xm,n of the set [M+l, N] E Z, 
we denote by the same symbol the projector on the sites \M + 1, N], corresponding to the 
multiplication operator by xm,n- Therefore 

TrM,N _ TjM,N _ TT M,N _ j T M,N^. (o a\ 

v — XM,NU = U XM,N — XM,N<J XM,N- [o.b) 

We now consider two measures related to finite matrices as follows. 
Definitions: The measures dkM,N an d dkM,N on T ar & defined by 

f(e lX )dk M , N (X) := tr (f(V M ' N ))/(N - M) (3.7) 
f(e lX )dk M ,NW ■= tr (XM,Nf(U) X M,N)/(N - M), (3.8) 
for any continuous function f : S 1 — > C. 

Note that dkM,N is nothing but the counting measure on T associated with the spectrum of 
the finite block V M,N , and dkM,N that associated with the projection of U on [M + 1,N]. 
This former operator is unitary whereas the latter is not. 

We denote the trace norm by || • ||i and first show a slight generalization of |GTj allowing 
to get 

Lemma 3.1 With the above notations, assume 

\\(U M ' N - U) X m,n\\i = o{N -N), as N - M -> oo, (3.9) 

then 

N l T-,oo F^M (^'^ - tr (.XM,Nf(U)XM,N)) = 0. (3.10) 
Remark: 

The hypothesis is satisfied in particular if Ka,nk(U M,N — U) < oo and uniformly bounded 
in (N,M), as is the case with the definitions of U M,N above by means of (|3.3( I3.5j) 
Proof: 

We first note that it is enough to consider functions which are polynomials in z and z, z £ 
S\ Any / e CiS 1 ) can be approximated by trigonometric polynomials Pr — X^j=— r 
in such a way that if e > is given, there exists R(e) < oo so that 

sup -P fl(e) (0)|<e. (3.11) 

Hence we get using (|3.6|) . 

tr (f(V M ' N ) - XM,Nf{U)XM,N) = tr (xM,N(f(U M > N ) - f(U)) X M, N ) = 
tr (XM,N(P R (e)(U M ' N ) - P R (e)(U))XM,N) + 

tr (XM,N((f - PR(e))(U M ' N ) ~(f- P R (e))(U))XM, N ), (3.12) 

where the trace norm of the last term is bounded by 2e(N — M), so that it becomes 
negligeable when divided by (N — M). We are thus to consider z s and z s , with s E N. We 
can write for any s > 1 

s-l 
j=0 
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so that 

8-1 

Xm,n(U s ~ (U N > M ) S )XM,N = Y,WU\U - U n > m )xm,n(U NM Y~ j - 1 - (3.14) 

3=0 

Therefore, 

tr ( X m,n(U s ~ (U N > M ) S )XM,N) < s\\(U-U N ' M )x M ' N \\^ (3 _ lg) 



N — M - N -M 

The same result is true if s < 0, with all unitaries replaced by their adjoints. Thus, 
— R(e) < s < -R(e) and the hypothesis on the trace norm of (U — U ' )x M ' yield the 
result. M 
Then, restoring the dependence on u in the notation, we get by the same arguments as 
in the self adjoint case, that the density of states is almost surely the limit in the vague 
sense of the measures dkM,N and dkM,N as ./V — M — > oo. A proof is provided in Appendix 
for completeness. 

Proposition 3.1 For any continuous function f : S 1 — > C, 



N— M— >oa 



T JT 



lim / f(e l *)dk M , N (\) = / f(e l *)dk(\) a.s. , (3.16) 



and the support of the density of states dk coincides with S, the a.s. spectrum ofU^ 



4 Thouless Formula 

The link between the density of states and the Lyapunov exponent is provided by an 
analysis of the spectrum of the finite unitary matrices V M)N . It reads 

Theorem 4.1 [Thouless Formula ] For any z £ C \ {0} 

j(z) = 2 [ ln\z-e iX '\dk(X') +ln(l/t 2 ) -ln\z\. (4.1) 
Jt 

Remarks: 

0) The identity ^(1/z) = ^(z) holds. 

1) It follows from the above formula, as in Theorem 4.6 in |(tT| . that the integrated density 
of states is continuous and satisfies 



\N{\ X )-N{M)\<^^^^ V where N{\) = jjk{\% 

by an argument of Craig and Simon |CSj . 

ii) In case z = e lX G S l , the formula can be cast into the form 



(4.2) 



7 (e jA ) = / ln(sin 2 ((A - A')/2))dA:(A / ) + ln(4/t 2 ), (4.3) 
from which we recover the estimate < 7(e* A ) < ln(4/t 2 ) that follows from the form of the 



transfer matrices 1)2.10)1 . 

The proof of this version of Thouless formula is given at the end of the section and 
we proceed with a Corollary and an application of this formula. The Corollary essentially 
expresses the radial derivative of the Lyapunov exponent as the Poisson integral of the 
density of states measure dk, which allows to recover the a.c. component of dk by a 
limiting procedure. 
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Corollary 4.1 For any e > and any X' G T, 

lim 7(e iA V e ) =7(e iV ), (4.4) 

| 7 (e iA 'e ±e ) = T jf |i^$!^*(A) - T^K^)- (4.5) 

Therefore, if n(X)dX/2ir denotes the a.c. component ofdk(X), 

lim |"7(e a, e-) = n(A') = ^l(e a '), (4.6) 

where the limit and the derivative exist for Lebesgue almost all X' G T. 
Remark: 

As in |CS| . it follows also from the subharmonicity of 7(2), that if j(e tX °) = 0, then 

7 : S 1 — > M + is continuous at e* A °. 

Proof: 

Let us first consider the second statement with lower indices only. We compute 

7 (e iX 'e- e ) = e + ln(l/t 2 ) + f ln(l + e~ 2e - e~ e 2 cos(A - X'))dk(X), (4.7) 

Jt 

which we can differentiate under the integral sign as long as e > to get 

= 1 + { -"^■" ! "'"!#) 

de n ' 7 T 1 + e _2<E — e~ e 2 cos(A — A ) v; 



1 — p~ 2e 

-dfc(A) =P[dfe](e iA 'e- £ ). (4.8) 



1 + e~ 2e - e" £ 2cos(A - A 

The existence for almost all A' G T of the limit and the first equality in (|4.6|) is a direct 
consequence of the above equality. The existence and equality with the derivative at zero 
for such A' follows from the mean value Theorem. To get the first statement, notice that 
1 + e~ 2e — e -e 2cos(x) > 2e _e (l — cos(x)) in formula (|4.7jl above yields 

0< -ln((l + e- 2e -e^2cos(A- A'))/4) < - ln(2 e - e (l - cos(A - A'))/4) = 
e-ln((l-cos(A- A'))/2), (4.9) 

where the last function is in L 1 (T,<i/c) by Thouless formula. Therefore, an application of 
the dominated convergence Theorem shows we can take the limit e — > inside the integral 
to get the result. I 

We consider now the properties of characterized by i.i.d. phases 9% and a£ in the 
definition 1)2. 2|) . assuming one set of phases is uniformly distributed on T. In that situation, 
not only can we can prove the transfer matrices have a (positive) Lyapunov behaviour, 
but we can also exactly compute the Lyapunov exponent 7(e lA ). This shows that in this 
situation, the spectrum of is almost surely singular, in view of the unitary version of 
the Ishii-Pastur Theorem proven in [BHJj. This strengthens the corresponding results of 
BHJ , Theorem 4.1 and Propositions 5.4. There Furstenberg's Theorem is applied to prove 
positivity of the Lyapunov exponent, so that no value for 7 (e* A ) is provided. 
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Theorem 4.2 Let (O^)kez an d {ct%)k& be i.i.d. on T and assume the distribution of either 
the 9% 's or the ott 's is uniform on T. Then, for any A 6 T, 

dk{\) = d\/2iT, and 7 (e a ) = ln(l/t 2 ) > 0, (4.10) 

therefore, 

o'iU^a.c = and o-(Uu) s i ng . = S 1 almost surely. (4-11) 
Remark: 

The assumption on the distribution of the phases actually implies that the ry^s are i.i.d. 
and uniform on T, see Lemma 14 . 1 1 b elow . This explains why the a.s. spectrum coincides 
with S 1 and why the density of states is flat. 
Proof of Theorem 14.21 

We first use the following lemma of purely probabilistic nature proven in Appendix. 



Lemma 4.1 Under the hypotheses of Theorem \4-^ the 's are i.i.d. and uniform on T. 

Then we show the density of states is uniform for uniformly distributed phases. Expanding 
()2,2|) of the r/fc(o;)'s we can write for any 

k=ki,k2,---,k n -i 

J^exp -i^Pimiu) (U ) jM (U ) klM • • • (t/o) fcn _i,j, (4.12) 
1 \ fe£ / 

where Uq corresponds to when all phases % = and where £ is a finite set of indices 
depending on j,k,n and pi are integers. Observing that the variables ^(u^'s all appear 
with the same sign in ()2.1|) . no compensation can take place between contributions of 
different matrix elements above and one at least among the integers pi, for / £ C is stricly 
positive when Using independence and the characterization K(e~ imVk ^') = 5 m fi of 

the uniform distribution, we get 

= / e inX dk(X) = 6 n , (4.13) 

Jt 

and the first statement follows. The second equality is a consequence of Thouless formula 
together with the identity 

In |1 - e iX \d\ = 0. (4.14) 



o 

The singular nature of the almost sure spectrum of U u comes from the unitary version 
of Ishii-Pastur Theorem proven as Theorem 5.3 in BHJ , which is independent of the 
properties of the common distributions of the ct^'s and #fc's and only requires ergodicity. 
Finally, Pr op osition 13 . 1 1 yields the result about the support of the a.s. singular spectrum. 
■ 

We compute here, for the sake of completeness, the density of states and Lyapunov 
exponent for the deterministic free operator Uo corresponding to XJ W in case ry~ = 0, VA; 6 Z. 
In this case, equation (|3.16|) of Proposition 13.11 becomes a definition of the free density of 
states dko, provided the limit exists. That the limit exists, is the content of the next 
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Lemma 4.2 The free density of states dko exists when defined for any f £ C(S' 1 ) by 

[ f(e lX )dk (X) = lira / f(e iX )dk M , N (X). (4.15) 

As we know essentially everything about the purely a.c. operator Uq, we can also use 
a direct approach to perform these computations. In particular, the integrated density of 
states of Uo can be defined as the distribution function on T of the band functions yielding 
the spectrum So of Uq. This direct approach of the density of states coincides with the 
above definition, see the proofs of Proposition 14. II and Lemma 14.21 in Appendix. We note 
here that the spectrum of Uo consists in the set 

S = { e ±*( arccos (r 2 -* 2 cos (y))) ; y g t}. (4.16) 

We get in particular that So is the support of the density of states whereas S§ is that of 
the Lyapunov exponent: 

Proposition 4.1 If No, dko and 70 denote the integrated density of states, the density of 
states and Lyapunov exponents of Uo, respectively. We have for A £ T ~] — 71", it], 

f 1 J si " (A)l w dX if |A| < arccos(r 2 - t 2 ) 

dko{\) = { 27rV*Mr 2 -co S (A))2 11 K ' (4.17) 

y otherwise 

_L arccos C r2 -™W ) if A € [- arccos(r 2 - t 2 ), 0] 



NoW = < V / 2 , ' (4.18) 

h-i arccos ( r ~T (A M if A G [0, arccos(r 2 - t 2 )} 



( if |A| < arccos(r 2 — t 2 ) 

7o(e 4A ) = , -i (r*-co s{ x)\ (4-19) 



^ cosh f 12 J j otherwise. 

Finally, Thouless formula \4-l\) holds true for these quantities with z = e tX , A G T. 
Remarks: 

Note that the density of dko(X) diverges as l/y|A — arccos(r 2 — t 2 )| at the band edges and 
behaves as l/2irt as A — > 0. 

The integrated density of states iVo(A) tends to its values and 1 as y/\X — arccos(r 2 — t 2 )\ 
at the band edges. 

Also, in keeping with the fact that Uo becomes a shift if t = 1 and the identity as r = 1, 
iVo(A) becomes linear in A as t — ► 1 and a step function as r — ► 1. 
The Lyapunov exponent, where non zero, is equivalently given by 



, i r 2 — cos(A) \ ( r 2 — cos(A) \ 2 , 
70 (e iA )=ln[ TT^ + Vl "M- M20 ' 



It is an even C°° function of A on {|A| > arccos (r — t )}, strictly increasing on [arccos (r — 
t 2 ),7r]. And djo{e lX )/dX behaves as 1/yA — arccos(r 2 — i 2 ) as A — > arccos(r 2 — i 2 ) + . 
Given Lemma 14.21 above, it is clear that Thouless formula holds for the above quantities. 
A direct proof of this fact is nevertheless given in Appendix. 
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4.1 Proof of Thouless Formula 

We now turn to the proof of Theorem 14.11 Writing down explicitely the effect of the 
boundary conditions at N > M on the coefficients of the eigenvector 1)2. 7 Jl we obtain the 
following relations, which depend on the parity of N and M. Let ip M > N = XM N^P and 
consider 

V M,N^M,N = e iX rt pM,N in /2 [M+1)A ^ (421) 
We get by inspection, 

Lemma 4.3 Assume j4.21\) is satisfied. Then, if M is even 

c M+2\_ r , h( i\^_ n 1 / -it(r-e~ tX ) \ , . 

cm+3 ) t 2 \(r - e lA ) + r(r - e lA ) ) 

If M is odd, 

rw+J ^ - c M+ ib 2 (e iX ) = c M+ i\ ( a ) ■ (4.23) 



. CA/+2 ) it 

Similarly, if N is even, 

Civ-i J " CNb3(e ) = CN ¥{ -it(r - e~ iX ) ) ■ (424) 
If N is odd, 

These relations together with the formulas (|2.12|) allow to describe the spectrum of V M,N 
in a convenient manner. 

Corollary 4.2 Lei M < N be fixed and consider non zero vectors a\,a 2 € C 2 smc/i i/iat 
aj (e iX ) e (b j+2 (e iX )C) L , j = 1,2. Then, e iX e o"(F M ' JV ) ijf 

(ai(e iA )|T(A^/2-l)---T(M/2 + 2)6i(e iA )) = 0, M, TV even 

(a 2 (e iX )\T((N + l)/2 - 1) • • • T(M/2 + 2)&i(e iA )) =0, M euen , iV odd 
(ai(e iA )|T(A^/2-l)---r((M + l)/2 + l)o 2 (e iA )) = 0, M odd ,N even 
{a 2 (e iX )\T((N + l)/2 - 1) • • • T((M + l)/2 + l)o 2 (e iA )) =0, M, iV odd (4.26) 

Remark: 

In particular, a possible choice for the a,-' is 

ai (e lA ) = 6!(e- iA ), a 2 (e iA ) = 6 2 (e~ lA ). (4.27) 

Each of the above quantities denotes a matrix element of a product of transfer matrices of 
the type ()2.12|) . which depend on e lX , and will be linked in the limit N — M — > 00 to the 
Lyapunov exponent. 

Let e %x = z 6 C \ {0} and no, TUq £ Z. Defining 

= T(n - 1) • • • T(m + 2), (4.28) 

one sees that the matrix elements {aj{z)\& mo,no (z)bk{z)) correspond to those in the above 
corollary for values N = 2n$,N = 2uq — 1,M = 2mo,M = 2rriQ + 1, depending on the 
choice of indices j,k. 
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Lemma 4.4 For any z G C \ S 1 and any indices j,k = 1,2 

1 



lim 



no-mo— too 2(uq — mo) 



ln\( aj (z)\^ m ^(z)b k (z))\ 



ln\z-e iX '\dk(X') + ln(l/t) - ln(|z 



1/2 x 



(4.29) 



Proof: We note that for any k G Z, there exist 2x2 matrices A(k) , B (k) , C (k) such that 
(with 2 = e iA ) 



T(ife) = zA(fc) + B(k) + C(Jfe)/z, where A(k) 







(4.30) 



(k) (k) 

Also, for any j = 1, 2, there exist vectors bj , a - , k = —1,0, 1 such that 

( \ (i) , (o) , (-i) / 

a k (z) = za k ' + a k ' + a k J /z, 

b k (z) = zb^ + bf + b k ~ 1) /z, 



(4.31) 



M) - „(1) 



are the only non zero vectors with the choice (|4.27[) . Thus, taking 



where b 

into account the above property , 

P hk {z) = z^^+^(a j {z)\^ m °^{z)b k {z)) 



(4.32) 



is a polynomial in z of degree 2(uq — mo) + 2 — (A: + _?'). Let Pj tk be the coefficient of the 
highest power of z of Pj tk . Then, because of corollary I4.2[ we can write 



deg P jtk 

P j>k (z)=p j>k ]"] (z-e lX '), 



(4.33) 



1=0 



where {e lXl } is the set of eigenvalues of V M,N and we compute 



M.N 



1^1 = 1(^1 



^2(no-mo) 



n A(i)bF)\= 

l=mo+2 

o\ ( n °~ m °)~ 2 

lj 



-it 



^2(no-mo) 



(4.34) 



where Kq, K\ are some constants that depends on j, k and t. Therefore, for any z £ C\S 1 , 

deg P jik 



no-mo— »oo (no — mo) 



lim 



ln(l/t 2 ) + lim V 

rt.r\ — mn — ^nn ' ^ 



In z — e 



/.':)-;• ii— x ^ (no — mo) 



(4.35) 



Introducing the continuous function f z : S 1 —> R given by /^(x) = In \z — x\, the last term 
can be written 



deg p. 



//I,-//',)— 'X ^— - no — mo M-N— >oc N — M 



by application of Lemma l3.1l and Proposition l3.11 This ends the proof of the lemma. 
Then we make use the following easy lemma 
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Lemma 4.5 If <I> : C 2 — > C 2 is linear and aj, bj € C 2 , j = 1,2 are such that span (ai, 02) = 
span (61,62) = C 2 , t/ien [|$[| := max^fc |(aj |<&6fc) | is a norm for $, 

noting that its hypothesis is satisfied by a^(z), bj(z), for all z 7^ —1, and of the fact that 
the Lyapunov exponent is defined independently of the norm used in (|2.13|) to deduce 
that (|4.29|) actually equals half the Lyapunov exponent. Finally, the fact that both the 
Lyapunov exponent and the right hand side of l|4.29j) are subharmonic and coincide on 
C \ S 1 implies the relation (|4.1|) on C as well, by classical arguments, see |CS| . This ends 
the proof of Thouless formula. I 



5 Properties of the Density of States 

We mentionned several times the analogy between our unitary operator U w and Jacobi 
matrices corresponding to the self-adjoint case. In this section we slightly drift away from 
the physical motivations underlying the study of (|2.1|) and consider more closely the links 
between these cases. The analogy is made clearer by the following Lemma which will be 
useful later. 



Lemma 5.1 Denoting unitary equivalence by ~, we have 
Uu ~ D^Sq, with D u = diag {e' 1 ^ } 



(5.1) 



and 



I 



S'n 



rt 
rt 



-t l 
-rt 

-tr 



\ 



V 



rt 
r 2 
rt 

-t 2 



-rt 
-tr 



Uq, 



(5.2) 



/ 



-i 2 , k £ 



where the translation along the diagonal is fixed by (^2fc-2|5'o ( / J 2fc) 
Remarks: 

In some sense, the Lemma says that, up to unitary equivalence, U w is a unitary analog of 
the one dimensional discrete random Schrodinger operator where the a.c. unitary Sq plays 
the role of the discrete laplacian, the pure point diagonal operator plays the role of the 
potential on the sites, and the operator sum is replaced by a product. 
We also recall that tridiagonal unitary matrices are spectrally uninteresting as they either 
correspond to a shift of to infinite direct sums of blocks of size one or two, see Lemma 3.1 
in [BIO] . 

The Lemma also shows that our operator \J W is essentially a product of an absolutely 
continuous unitary and a pure point unitary, whereas it was constructed in Section 2 of 
|BH.Ij as a product of two pure point unitaries. 
Proof: 

Let us define a collection of rank two operators by 

p j = \<Pj){<Pj\ + IVj+iXVj+il. JGZ, (5.3) 

and the unitary V by the direct sum 

V-e / ir t 
1. P2 ^[- it r 



V 



2.7-1- 



(5.4) 
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It is just a matter of computation to check that we can write 

U u = (UuU^Uq = V- l D u VXJ Q = V^DUVUoV-^V = V^ 1 (D aJ S )V, (5.5) 
with the required properties for Sq and D u . M 

Now, forgetting that the phases r/^ are in general correlated random variables, see (|2.2|) . 
if we consider them as i.i.d., but not necessarily uniformly distributed on T, we get some 
unitary Anderson-like model. This is where we depart from the physical motivation, as it is 
recalled in Lemma 4.2 in |BHJj that independence of the rj^'s is associated with a uniform 
distribution. 

5.1 Support of the Density of States 

Nevertheless, assuming the random phases {rf£}k& are i-i-d- according to the measure d[i 
on T, we can characterize the almost sure spectrum of U u in term of the support of \x and 
of the spectrum So of Uq. 

Theorem 5.1 Under the above hypotheses, the almost sure spectrum ofU^ consists in the 
set 

£ := exp(z supp//)So = {e ia T,Q \ a G supp//}. (5-6) 
Remark: 

In the case where the %(w) are i.i.d. and uniform on T, we recover the fact that the almost 

sure spectrum of U w is S 1 . 

Proof: 

To show that E belongs to the almost sure spectrum, we simply construct Weyl sequences 
corresponding to the corresponding quasi-energies, with probability one. We know from 
Section 6 of |BHJj that for any e lX G So, there exists a generalized eigenvector ip\ such 

= J2 c i( A )^'' = e ^A, and < K < \ Cj (\)\ < 1/K, Vj G Z, (5.7) 

for some K > 0. The last property can be checked also by means of the transfer matrices 

Let a G supp/x. Then, for all e > 0, there exists a set I t B a such that |J e | < e, and 
n(I e ) > 0. With the notation u(k) = r]k(u>), k G Z, we define for all n G N and k G Z, 

A n (k) = {uj(kn) G I e , co(kn + 1) G I e , ■ ■ ■ , u;(kn + n - 1) G J e }. (5.8) 

Due to the assumed independence, we have for any k, F(A n (k)) = (i(I e ) n > so that for 
any n > 0, by Borel-Cantelli, P(U fcG z-4n(^)) = L 

Let A n (k) = {kn, kn + 1, • • • , kn + n — 1} denote the set of indices appearing in A n (k) 
and consider now 

VVc(A)= Yl c i( A )^' = x(An(fc)h/>(A), (5.9) 

jGA n (fc) 

where x(A n (k)) is the projector on the span of {^>j}jeA n (k) Because of (|5.7j) . 

Uo^kW = e a Vn,fc(A) + ffo(A) + fl+ , (5.10) 
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where the vectors i? • have at most four components close to the index j and 

W-R-fW — ^) where R is uniform in j. (5-H) 

Also, by construction of A n (k), Uq and U^, we have 

||^n,fc(A)-e iQ C/ ^(A)|| < WiU^ - e ia U )x(A n (k))\\\^ k (X)\\ 

= 0(e)||^ n , fc (A)||, (5.12) 

where the estimate 0(e) is uniform in n and k. Therefore, for all e > and all n > 0, there 
exists, with probability one, a k such that A n {k) and the corresponding ip n ,k(^) have the 
above properties so that 

\\U^ n>k (X) -e i ( Q+A )Vn,fc(A)||/||^, fe (A)|| = 

(||(^ - e ia C/ )V'n,fc(A) + e ia (U - e iA )Vy fc (A)||)/||V>n,fc(A)|| < 

0(e) + 2i2/||^ n>fc (A)|| = 0(e + 1/n). (5.13) 

It remains to chose n = [1/e] to conclude that ^i a +^> G cr(U u ) almost surely. 

Let us now show that S 1 \ £ belongs to the resolvent set of U u . In order to do so 
we use Lemma 15.11 Therefore, we can consider as well the spectrum of the product D^Sq 
to which the perturbation theory recalled in Chap.l, §11 of |Yaf| for example, applies. 
In particular, dropping the u in the notation as randomness plays no role here, if we 
know that for all j G Z, rjj G [a, (3] C T, then cr(D) C (61,62) where (61,62) denotes 
the corresponding arc on the unit circle swept in the positive direction from 61 G S 1 to 
62 G S . We denote by |(<5i,^2)| the length on the torus of this arc. Since cr(So) = ^0 
corresponds to the symmetric arc (e~ l arccos ( r _ * ) e «arccos(r -* )) ) perturbation theory tells 
us that after (multiplicative) perturbation by So, the spectrum of U ~ DSq is a subset of 
an arc of wider aperture than (^i,^)- Quantitatively, Theorem 8, p. 65 in Yaf tells us 
that the arc ( e * arccos ( r2 -* 2 )^ 2) e -«arccos(r 2 -i 2 )^^ b e i on g S to the resolvent set of U, provided 
|(<$l,<$2)| < I (e 1 arccos ( r -t ) ) e -«arccos(r -t rp^jg condition simply insures that the subset 
of the resolvent set we are talking about is not reduced to the empty set. This is enough 
to get the result in case the support of \x is such that £ is connected. In case this set 
is not connected, as |Eo| > 0, it consists of a finite set of connected components, each 
of which can be associated with the convex hull of sufficiently far apart subsets of the 
support of fj,. Denoting these subsets by rrij, j = 1, ■ ■ ■ , N and the associated arcs on 
S 1 by (M\(j) , M2(j)) , we have that the spectrum of D is the disjoint union of subsets 
cjj satisfying Oj C (M\(j), M^fJ))- The same argument as above says that the spectrum 
of DS is confined to the finite union of arcs ((e iarccos ( r2 -' 2 )Mi(j), ( e - iarccos ( r2 -* 2 )M 2 (j)), 
which ends the proof of the Theorem. I 



5.2 Analyticity of the density of states 

At the price of some combinatorics, we can further exploit the relation (|4.12|) in order 
to obtain a condition on the common distribution of the rj^s ensuring the analyticity of 
the density of states. Recall that a function / on T is analytic, if and only if its Fourier 
coefficients / satisfy an estimate of the form 

\f(n)\ < Ae~ Blnl , Vn G Z, (5.14) 

for some positive constants A,B. We have 
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Theorem 5.2 Assume the rjk's are distributed according to a law that has an analytic 
density f characterized by the estimate {5.1$ with A,B > 0. Then, if 

B > ln(l + 2rt) + In A, (5.15) 

the density of states dk admits an analytic density, so that the integrated density of states 
N is analytic as well. 

Remarks: 

Asf(0)=J T f(r 1 )dr ) = l,A>l. 

When the Theorem applies, it prevents the Lyapunov exponent from being zero on a set 
of positive measure. 

This result has to be compared with the Proposition VI. 3.1. of |CL| stating a similar 
result for the d-dimensional Anderson model. 



As an immediate consequence, using r 2 + 1 2 = 1, we get the following 

Corollary 5.1 If the % 's have an analytic density f, characterized by \5.14\) with B > 
ln^4 ; then there exist r + {f) and r~~(f) in ]0, 1[ such that the density of states is analytic 
provided the reflexion coefficient r satisfies 1 > r > r + (f) orO < r < r~~(f). If B > ln(2A), 
The density of state is analytic Vr 6 [0, 1]. 

Remark: 

It is easy to check that in both the extreme cases r = 1 and r = 0, the density of states is 
analytic. Indeed, if r = 1, dk(X) = f(X)dX, where / is the density of the n^s, whereas if 
r = 0, dk(X) = dX/(2n). 
Proof of Theorem 15.21 
By hypothesis, for any n 6 Z, 



e ir > n f(r))d V 



<Ae~ B ^. (5.16) 



Then, in (|4.12|) above, YliecPl = n i so that using independence 

mVi\U>i)\ <A n e- Bn Yl \(Uohk 1 \\(Uo)k 1 M\---\(Uo)k n . 1 j\ (5.17) 

fel,fc2,-",fcn.-l 

Here the sum carries over a set of indices that form paths of length n + 1 from index j 
to index j. The allowed paths are those giving rise to non zero matrix elements (t^o)/,m 
in the sum above. In order to compute this last sum, we proceed as follows. Let us 
introduce more general j-dependent subsets C n -\{j) of indices of Z™ -1 that appear in the 
computation of the matrix element ((po\U™(pj). This set consists of paths of the form 
{&o = 0, k\, fo, • • • , k n — i, k n = j} of length n + 1 in Z from to j with the condition that 

k m +i - k m € {0, +1, -1, +2} if k m is odd 

k m+1 - k m G {0, +1, -1, -2} if k m is even, (5.18) 
for all m = 0, 1, -, n — 1. Let us define 

Sn-i(j):= £ l(^o)o,fc 1 ||(^o)fc 1 , fc3 |---|(^o)fc B _ ll il, (5-19) 
C„_i(0) 

where the matrix elements |(£/o)zro| are given by r 2 ,rt and t 2 respectively, when \l — m\ 
equals 0, 1 and 2 respectively. This quantity actually gives a crude upper bound on the 
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probability to go from site to j in n time steps, under the free evolution. It is crude in 
the sense that it does not take the phases into account during that free evolution. 

We are actually interested in the computation of S n -i(0) and of the similar quantity 
appearing in the computation of (<pi\U™(pi), which correspond the the sum in the right 
hand side of ()5.17|) . in the asymptotic regime n — > oo. The case of the matrix element 
{<Pi\U™tpi} being similar, we only consider S n _i(0). 

The plan is to use a transfer matrix formalism to evaluate the generating function 
associated with S n -i(j) and then to compute the asymptotics of 5 n _i(0). In view of 
(|5.17j) . the following proposition implies the Theorem. 

Proposition 5.1 For some constant c > 0, 

5 w _i(0) = 1 ' (l + o(l)) asn^oo. (5.20) 
\/n 



Proof of Proposition I5.lt 

Let 

Pn{x)= S n-lU> j (5.21) 

-2n<j'<2n 

be this generating function which we split into two parts P n (x) = P^(x) + P~(x) where 
P±(x) = Y, S n-iU> j - (5-22) 

-2n < j < 2n 
even 



J 



odd 



Clearly we have for n = 0, 1, 

P+( x ) =r 2 ,P~(x) =0,P+(x) =r 2 +t 2 x~ 2 ,P^(x) =rt{x + x~ 1 ). (5.23) 
It is readily shown by induction that a transfer matrix allows to compute P n (x) for any n: 



Lemma 5.2 For any n > 0, 

( Pn + i (x)\ = (r 2 + t 2 x' 2 rt(x + x- 1 )\( P+{x) \ 
\P~ +1 (x) J \rt(x + x- 1 ) r 2 + t 2 x 2 J \P-(x) J 1 

with P + (x) = r 2 ,P "(x) = 0. 

Denoting by T{x) the transfer matrix defined in this Lemma, and introducing the parameter 
t = t/r G]0, oo[, (5.24) 
we rewrite it as 

T{x)=r 2 (\ +T2x ~l r 1 (l + r a l) V (5.25) 

We will consider first the case t ^ r <J=^ r / 1. The case r = 1, for which more can be 
said about S n -i(j), see Proposition 15.21 is dealt with below. 



17 



5.2.1 Case r ^ 1 

The eigenvalues of T(x) are given by r 2 times \±(x), where 

A±(x) = {l + r(x 2 + x" 2 )/2 ± v/(i + r(x 2 + x" 2 )/2) 2 - (1 - r 2 ) 2 } , (5.26) 
so that 

T"(x) = r 2 M(x)( A + (x) A ^^(x)- 1 (5.27) 

with 

j4f3 . ) = / A +(^)-( 1 + ^ 2 ) A_(x)-(l + r 2 x 2 )\ 
1 ' \ r(x + x _1 ) r(x + x _1 ) /' 1 ' J 

For the moment, x is just book keeping parameter, so that we ignore the potential problems 
of the definition of A(x) in case the eigenvalues are degenerate and we further compute 

?(«)) _I " (l) (o)" (5 ' 29) 



r 2n r(x + x 1 ) 



2^(1 + r(x 2 + x" 2 )/2) 2 - (1 - r 2 ) 2 

A+(x) n+1 - A_(x) n+1 - (A+(x) n - A_(x) n )(l + r 2 x 2 ) ' 
r(x + x- 1 )(A+(x) n - A_(x) n ) 

We note at this point that one checks, using the binomial Theorem, that despite the 
presence of square roots in the expressions for P^(x), these quantities actually are given 
by finite Laurent expansions in x, as they should. Focusing on P+(x) we can rewrite with 
the shorthand \f- for the square root of the denominator above 

P+(x) = (5.30) 
>J " Ti ' ' ' '"' 1 ({\ + {xY - \-(xT)4(x- 2 + x 2 ) + ^(A + (x)" + A_(x)») 



2^T- V 2 2 

The quantity of interest to us is S' n _i(0), the coefficient of x° in the expansion of P^{x). 
Substituting e td for x in -P^ - , we get a trigonometric polynomial whose zero'th Fourier 
coefficient is obtained by integration 

5„_i(0) = / P+(e w )d8/(2n). (5.31) 

It remains to perform the asymptotic analysis of the above integral as n — > oo. It is 
a matter of routine to verify the following propereties: The eigenvalues, as functions of 
6 G T ~] — 7r,7r], are continuous. If r < 1, they are real valued, with discontinuity of the 
derivative at 9 = ±ir/2, where they cross and are given by 1 — r 2 . At all other values of 9, 
they are C°° and they satisfy 

\+{e ie ) > A_(e ie ), with \+{e ie ) > 1 - r 2 . (5.32) 

If r > 1, the eigenvalues become complex conjugate. Let 9 C = arccos( T ^ 2 )/2 be the 
critical value where the square root becomes zero. If 9 G [9 C ^ — 9 C ] U [— tt + 9 C , —9 C ], the 
eigenvalues are complex conjugate, of modulus |1 — r 2 |. Otherwise they are real valued, 
and satisfy (|5,32|) as well. Therefore, the asymptotics as n — > oo of (|5,31|) is determined 
by A+ only. Moreover, in both cases, ln(A+(e i61 )) admits non degenerate maxima at 9 = 
and 7r, where A_|_ reaches its maximum value (1 + r 2 ). Therefore, Laplace's method yields 
the asymptotics of the Proposition. I 
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5.2.2 case r — 1 

The course of the proof being the same, it is presented in Appendix. However, instead of 
computing 5 n _i(0) as n — > oo, we can get exact forms for all 5 n _i(j)'s. The Proposition 
we actually show is 

Proposition 5.2 

1 / 2n — 1 \ 

Sn-i(j) = ^( y 2 + J, ~2n < j < 2(n - 1), j even 
1 / 2r? — 1 \ 

S„-i(j) = 2^ (j -_ 1)/2 + n J, -2n + l<j<2n-l, j odd (5.33) 
Remark: 

Of course, Stirling's formula for n large yields vrovosition \5.1\ with r = t = \j\pl: 
2 n \ n / Jnri 



6 Appendix 

Proof of Proposition I3.ll 

We have by definition, 

r 1 N 

f(e iX )dkl jN {\) = ——- Y, (<Pj\f(P»)<Pj), (6-1) 

Jj j=M+l 

where, depending on the parity of M and N and due to the fact that / is uniformly 
bounded, the right hand side can be rewritten as 



1 / m \ 1 

E (V2k\f{u^2k) + (^+11/(^)^+1) + Ofijj—jj) = 

\fc=(M+l)/2 / 
N/2 \ 

E ^o|/(^h)^o) + <^i|/(^H Vi> + O/f^). (6 

M+l)/2 / 



v fc=(M+l 

Now, by Birkhoff ergo die theorem, there exists of measure one such that for all u G £lf, 

iviooiV^M £ = ^E((^|/(^)^)),Vi G Z, (6.3) 

fc=(M+l)/2 



therefore, 



' tr (x MlA 7(^)) - ~ m<PoW u )<po) + ^o|/(^)^))) • (6.4) 



N — M v/v J v " 2 

Then, C(S' 1 ) being separable, we have the existence of a countable set of {fj}j£N, dense 
in C(S' 1 ), for which the above is true, on a set of probability one, which proves the almost 
sure convergence stated in the proposition. 

Now assume e lX ° S and take a continuous non negative / such that f(e lX °) = 1 
and /|s = 0. Then f(U u ) = a.s. so that / f{e iX )dk{\) = and e iX ° # supp k. 
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Conversely, if e iXo supp k, there exists a non negative continuous / with f{e lX °) = 1 and 
/ f(e lX )dk(X) = 0. Hence, a.s., (v?o|/(t^)<£o) + {v\\f{Uw)<P\) = 0, therefore, by ergodicity, 
(fjlfiU^tfj) = a.s. for any j and f{U w ) =0. As / is continuous and equals one at e lX °, 
we get that e iX ° £ E. ■ 

Proof of Lemma I4.lt 
We only deal with the case where the #£"s are i.i.d. and uniform, the other case beeing 
similar. Let $m(ra) = E(e mj? fc) be the characteristic function of the random variable rf£, 
and similarly for a£, and $g(n) = <5 nj o- Then, using independence, 

^(n) = $ (n) 2 $ a (n)$ a (-n) = 5 nfi \ <S> a (n) | 2 = <5„, , (6.5) 

so that the %'s are uniformly distributed. Consider now 

$v ko m v ~,v kj («o,ni, -, nj ) = ECe^U^). (6 . 6 ) 

We can assume the fej's are ordered and we observe that % and rjk+j are independent as 
soon as j > 2, see (|2.2j) . Therefore, we can consider consecutive indices ki and deal with 

^MH.r.l+it' 1 !.^,-,^) = (6-7) 
]W e ino0fc_i+i(no+ni)0 fc H Mfaj-j+nj^fc+^-i+njflj)^ jf a 

where the second expectation contains a^s only. Then 

<S> e (n )<S> e (n + m) ' ' ' + «j) $ e(%) E (/(a)) = 

<Wo<*ru,o • • • ^,oE(/(a, n)) = <^ ff E(/(a, 0)) = 6^, (6.8) 
whith the obvious notation, which yields the result. M 

Proof of Proposition I4.lt 

We first prove this Proposition with the definition of the density of states as the distribu- 
tion function of the "band functions" of Uq, to be defined below. Then we'll see in the 
course of the proof of Lemma 14.21 below the equivalence with the definiton as an average 
counting measure. The proof of Proposition 6.2 in |BHJj shows that Uq on Z 2 (Z) is unitarily 
equivalent to the operator multiplication by the matrix 



by the unitary mapping that sends i— ► e lkx /\^27r, and where L±(T) is the susbspace 
generated by even/odd harmonics {e lkx }k£i,- The eigenvalues of V(x) are 

\±(x) = e ±ia{x \ where a{x) = arccos(r 2 - t 2 cos(2x)). (6.10) 

We note that \±(x) = \±(—x) and 

V(x) = JV{—x)J where J=(? \ J . (6.11) 



1 0^ 

Hence, the corresponding eigenvectors x±(^) satisfy 

^(aOx±(z) = A±(a;)x±(aO and V(x)Jx±(-x) = A ± (x)Jx±(-a:), (6.12) 
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so that x±( x ) an d Jx±{~ x ) are linearly dependent. This is in keeping with the fact that 
the subspace of generalized eigenvectors is of dimension 2, see (|2.8[) . Also, one checks that 
for any phase /3 G] — arccos(r 2 — t 2 ), 0[U]0, arccos(r 2 — t 2 )[, 



= {x 1 ,x 2 , -x 2 - xi} C] - 7T,7r[. (6.13) 



Therefore, due to (|6.12j) . only half these points contribute for the computation of the 
density of states. We can now compute the integrated density of states No((3) as follows: 
Taking into account the normalisation by a factor 1/2-7T in the definition 1)3. lj) . the fact 
that supp k C [— arccos(r 2 — t 2 ), arccos(r 2 — t 2 )] and the symmetries, we have for any 
(3 G [-arccos(r 2 - t 2 ),0] 

if i r /2 

N o(P) = I d\X{-a(\)<(3<0} = ^ / ^X{cos(2A)>(r 2 -cos(/3))/t 2 } (6.14) 

x ,arccos((r- 2 -co S (/?))/t 2 ) j /r 2 - COs(/3) \ 

■ arccos I ^ . (6.15) 



2vr 7 2vr \ t 2 

A similar computation for (3 G [0, arccos(r 2 — i 2 )) yields ()4.18j) . Therefore, dko is absolutely 
continuous w.r.t. Lebesgue and, for any |A| < arccos(r 2 — t 2 ), dko(X) = N'(X)dX, from 
which the result on the density of states follows. In order to obtain the Lyapunov exponent, 
it is enough to observe that the transfer matrices 1)2,8)) T, now independent of k, are of 
determinant one and trace equal to 2(r 2 — cos(A))/t 2 . Therefore, it is readily checked that 
when the eigenvalues t±(A) of T 



r± (A) = (r 2 - cos(A) ± yV = cos(A)) 2 - t 4 )/t 2 (6.16) 
are complex conjugates, i.e. when |A| < arccos(r 2 — t 2 ), they are of modulus one, whereas 



max{|r+|, |r_|} = (r 2 - cos(A) ± x/(r 2 - cos(A)) 2 - i 4 )/t 2 , (6.17) 

if |A| > arccos(r 2 — t 2 ). It remains to use definition (|2.13j) to get 7o(e lA ). In order to prove 
the last statement, we first rewrite the right hand side of Thouless formula with dko(X') 
above as 

i [ l H(x-yn dx + ln2 (fU8) 



2vr J_ x vT 



x 



2 



by means elementary manipulations, changing variables to x = (r 2 — cos(A'))/t 2 and intro- 
ducing y = (r 2 — cos(A))/t 2 G [— 1, (r 2 + l)/i 2 ]. Hence we are to show that ()6.18j) above 
equals if y < 1 and ln(y + \Jy 2 — 1) if y > 1. We first deal with the case y > 1. We can 
differentiate (|6.18|) with respect to y under the integral sign to get 

1 f 1 dx 1 f dz 

- / = — / , (6.19) 

71" 7-i yl - x 2 (y - x) 27r 7 C yi - z 2 (y - z) 

where C is a contour in the complex plane surrounding the segment [—1,1] in the positive 
direction which does not contain y in its interior. By deforming the contour to a circle 
centered at the origin and of radius R > large enough, we pick a residue at y. As the 
integral on the circle is of order 1 /R, we eventually get in the limit R — > oo 

d ( 1 f 1 ln((x-y) 2 ) 1 1 

1 1 ' -dx + ln2} = (6.20) 



dy [2ir J_ x yT^ 2 J ^fy 
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as expected. The limit as y —* 1 + of the Lyapunov exponent fixes the constant to 0. Now, 
if y G] — 1, 1[, we first convert 1)6.18)1 to a contour integral along a path similar to the one 
above with te following difference. As the In is multivalued, with a cut from y along the 
real axis towards — do, the contour is attached to the point —1. By assumption, y does 
not belong to the contour of integration, so that we can now differentiate with respect to y 
under the integral sign and thus get the same contour integral 1)6.19)1 as above to consider. 
However, by expanding the contour to infinity, we get to residue this time, so that 1)6.18)1 
is constant for y G] - 1, 1[. As it is known f |GR| . # 4.224, p. 526) that, 



1 f 1 In x 2 , 2 f 1 In x , 2 /* / , 
— / ; v ' dx = - / - v 1 dx = - / ]n(wa.(t))dt = -In 2, 6.21 

we have, by continuity, that the integral is equal to zero on [—1, 1]. I 
Proof of Lemma 14.21 

We use freely the notations above. Let us introduce the eigenprojectors P±(x) associated 
with \±(x) such that 

V(x) = P + (x)X + (x) + P_(a;)A_(aj). (6.22) 

These quantities are analytic in x, in a strip including the real axis. Let / G C(S 1 ) and let 
us compute by means of (jUSJ) and the definition of L±(T) 

tr (XM,N\f( U o)XM,N) = ^2 ( t Pj\f(. U o)'Pj) = 

M<j<N 



j even 

M<j<N 



E h I ((J)|(/(A + (x))P + (x) + /(A_(x))P-(^))(J 

even \ / \ . 

E ^/(( 1 )|(/(A + (x))P + (x) + /(A_(x)) J P-(^))(j))^- (6-23) 



j odd 

M<j<N 



The summand being independent of j and uniformly bounded, we can rewrite the above 
trace as N — M gets large as 



f /(A+(x)) tr P+(x) + /(A_(x)) tr P_(x)dx + 0(1) = 
Jt 

f(X + (x)) + f(X.(x))dx + 0(l). (6.24) 



7V-M 

47T 

N-M 

Hence, with X±(x) = e^ %a ^ as in ()6.1U() . and taking into account the properties of a, we 
get 



Jt 



f(e iX )dk (X) = i- / /( e -(-)) + /( e --W)dr 



f7r/2 

| /( e M W) + /( e - M W)(ii, (6.25) 

°" -tt/2 



2tt 



which is easily seen to coincide with the " direct" definition of dko in the above proof. H 
Proof of Proposition 15.21 

As in that case a commun term ^ can be factorized, see 1)5.17)1 . we compute the generating 
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function of \C n -i(j)\, the cardinal of the set of relevant indices. Using the same symbols 
as above, we consider this time 

P n (x)= Yl \ C n-lU)W, (6-26) 

-2n<j<2n 

which we split into two parts P n (x) = P^(x) + P~{x) that satisfy for n = 0, 1, 

p+( x ) = l 7 P-(x) =0,P+(x) = 1 + x~ 2 ,p- (x) = x + x~ 1 . (6.27) 
As above, 

Lemma 6.1 For any n > 0, 

P n + +1 (x)\ = /l+X- 2 X + X- 1 } fP+(x) 
Pn+1 W Kx + X- 1 1 + X 2 J V Pn (*) 

with P + (x) = 1, P Q (x) = 0. 

By diagonalization of the corresponding transfer matrix, we get 

T n (x) = A(x)(° n ^W)"* 1 ( 6 - 28 ) 



[x + x 
where 

and we compute 

Using the binomial Theorem we obtain for P^(x) 

- x> a ( 2 r + -;) 

l=—n 

Pn(x) = E^ft+n 1 )' (6 ' 31) 

l=—n 

hence the end result. S 
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